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The interfacial instability and subsequent dynamics in a phase- separated two-component Bose- 
Einstein condensate with rotation symmetry are studied. When the interatomic interaction or the 
trap frequency is changed, the Ray leigh- Taylor instability breaks the rotation symmetry of the 
interface, which is subsequently deformed into nonlinear patterns including mushroom shapes. 
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I. INTRODUCTION 

The Rayleigh-Taylor instability [H-H (RTI) is an insta- 
bility of an interface between two fluids in a metastable 
state. For instance, when a layer of a heavier fluid is laid 
on a lighter fluid, the system is energetically unfavorable 
and the two fluids tend to exchange their positions. How- 
ever, if the two fluids are immiscible and their interface 
is flat, the exchange cannot occur without breaking the 
translation symmetry of the interface. Once an infinites- 
imal modulation arises on the interface, it exponentially 
grows due to the RTI, and the interface develops into 
complicated patterns such as a mushroom-shaped pat- 
tern. This kind of phenomena is found through nature 
in a wide scale ranging from laboratory to astronomical 
scales [5|. 

Recently, the RTI was predicted to be observed in a 
two-component Bose-Einstein condensate (BEC) [H, [?}• 
In Ref. [6[, a two-component BEC trapped in a tight 
pancake-shaped trap is considered, in which the two com- 
ponents separate into two semicircular shapes in the trap. 
When an external force is applied to each component in 
the direction to the other component, the RTI sets in 
and the interface deforms into a mushroom-shaped pat- 
tern. In Ref. [7], the initial state of a two-component 
BEC forms a domain structure in the axial direction of 
a cigar-shaped trap. The interaction between atoms in 
one component is then increased, and the RTI arises at 
the domain wall. 

The RTI is a symmetry breaking phenomenon, that is, 
even when the interface has symmetry (e.g. the transla- 
tion symmetry of a flat interface) an infinitesimal modu- 
lation grows exponentially and the symmetry is sponta- 
neously broken. However, in the above mentioned studies 
in Refs. [HQ, the symmetry breaking cannot be observed 
explicitly, since the relevant symmetry is broken from the 
initial state: the interfaces have edges or rims and have no 
translation symmetry. The inhomogeneous interfaces of 
the initial state strongly affects the deformation dynam- 
ics of the interfaces, obscuring the symmetry breaking 
nature of the RTI. 

In this paper, we propose a system to observe the sym- 



metry breaking RTI for a trapped two-component BEC. 
We consider a two-component BEC with rotation sym- 
metry, in which the two components separate radially 
and a "bubble" of inner component is surrounded by a 
shell of outer component. The interface between the two 
components has a spherical shape for an isotropic trap 
and a circular shape for a quasi-2D axisymmetric trap. 
If we change a parameter in such a way that the inner 
component tends to go out of the outer shell component, 
the RTI breaks the rotation symmetry of the interface 
and the spherical or circular interface is deformed into 
various patterns. The symmetry breaking RTI can thus 
be realized in a trapped BEC. The RTI that breaks rota- 
tion symmetry occurs in a variety of systems, such as 
supernova explosio n H , imploding targets in inertial- 
confinement fusion [10j, and collapsing cavitation bub- 
bles [ll|, [H. 

This paper is organized as follows. Section ITT1 provides 
a formulation of the problem and Sec. IIIII shows numer- 
ical results. Section IIII Al demonstrates the symmetry- 
breaking RTI and subsequent dynamics for an axisym- 
metric oblate system. Section IIII Bl shows dynamics for 
an isotropic trap and performs Bogoliubov analysis. Sec- 
tion IIII CI examines oscillation of interaction. Section IIVI 
gives conclusions to this study. 



II. FORMULATION OF THE PROBLEM 

We consider a mixture of two kinds of bosonic atoms 
with mass m\ and rri2 confined in trapping potential V\ 
and V2, respectively. The Hamiltonian for the system is 
given by 
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where ifjj is the bosonic field operator for component j 
and Ujjt is the interaction between atoms of components 
j and j' . In the mean-field theory, we assume that the 
atoms in each component occupy the same wave function 
ijjj and that the interaction potential is reduced to the 
Fermi pseudopotential, 

U jjf (r - r f ) = 27i?t 2 a jr (mj 1 + mj, 1 ^ S(r - r f ) 

= g jr 5(r-r f ), (2) 

where ajjt is the s-wave scattering length between the 
atoms in components j and j' . The system is thus 
described by the two-component Gross-Pitaevskii (GP) 
equation (j ± f), 
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The macroscopic wave functions are normalized as 
f\i/jj\ 2 dr = Nj with Nj being the number of atoms 
in component j. The two components are miscible for 
#n#22 > 012 and immiscible for gng 2 2 < g 2 2 - 

We solve the 3D GP equation (|3]) numerically using 
the pseudospectral method The initial state is the 
ground state prepared by the imaginary-time propaga- 
tion method, in which i on the left-hand side of Eq. (|3j) 
is replaced by —1. We then add a small noise to the ini- 
tial state as a seed that triggers the RTI. The dynamics 
do not depend on the initial noise qualitatively. 

In the following calculations, we assume a dual-species 
BEC with 85 Rb and 87 Rb, where the \f = 2,m f = -2) 
state of 85 Rb is component 1 and the \f = l,ra/ = —1) 
state of 87 Rb is component 2. This system has been re- 
alized by the JILA group [13], in which controlled phase 
separation was observed by changing the s- wave scatter- 
ing length an of 85 Rb using a magnetic-field Feshbach 
resonance, which is variable in the range an = 50-900aB 
with an being the Bohr radius. Since a 22 — 99<2b and 
ai2 = 213ae, the condition for the phase separation is 
satisfied for an < 458<2b- 
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FIG. 1: (Color online) (a)-(d) Dynamics of the col- 
umn density profiles D\ — f\ipi\ 2 dz (upper panels) and 
D2 = f \ip2\ 2 dz (lower panels) in an axisymmetric trap with 
{oj±,uj z ) — 2tt x (25,1250) Hz. The scattering length an is 
linearly increased from 80<2b to 240<2b between t = and 
t — 40 ms, and after that an is fixed to 240aB- The num- 
bers of atoms are N\ — N2 = 10 5 . The unit of the col- 
umn density is 10 12 cm -2 , (e) Cross-sectional phase profile 
<f>j — aig[il)j(z = 0)] of the lower half region of (c). The circles 
in (e) indicate examples of quantized vortices created under 
the caps of the mushrooms. The field of view is 65.4 x 65.4 
/im in (a)-(d) and 65.4 x 32.7 /im in (e). 



III. NUMERICAL RESULTS 

A. Rayleigh- Taylor instability in axisymmetric 
oblate systems 

We first demonstrate the dynamics for an axisymmet- 
ric oblate trap, Vj = mj\uj\{x 2 + y 2 ) + u 2 z 2 ]/2, where 
uj z ^> uj±. We assume that the gravitational sag is com- 
pensated and the two components share a common trap 
center. Figure [T] shows the time evolution of the density 
and phase profiles of the system, obtained by solving the 
3D GP equation ([3]) . The initial state is the ground state 
for an = 80<2b and Ni = AT 2 , which has the axisymmet- 
ric circular interface between the two components [Fig. [T] 
(a)]. The repulsive interaction of component 1 (inner) 
is then gradually increased and when it exceeds that of 



component 2 (outer) the system becomes metastable, i.e., 
the state in which component 1 surrounds component 2 
becomes energetically favorable. At t ^ 80 ms, the ax- 
isymmetry of the system is broken and the interface is 
modulated due to the RTI [Fig. (U (b)]. The modulation 
of the interface subsequently grows to become a four- 
fold mushroom shape [Fig. [1] (c)]. Quantized vortices 
are generated under the caps of the mushrooms in both 
components [circles in Fig. [1] (e)]. When the tops of the 
mushrooms reach the edge or the center of the system, 
a highly nonlinear pattern is observed [Fig. [1] (d)]. The 
n-fold mushroom shapes with n ^ 4 are also observed, 
where n is larger for a larger final value of an. 

The dynamics also depends on the ratio between the 
numbers of atoms N2/N1. Figure [2] shows the dynamics 
for N2/N1 = 1/9. After the repulsive interaction of com- 
ponent 1 is increased, the RTI causes modulation at the 
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FIG. 2: (Color online) Dynamics of the column density pro- 
files Di and D 2 for Ni = 1.8 X 10 5 and iV 2 = 2 x 10 4 . Other 
parameters are the same as those in Fig. [T] 



interface [Fig. |2] (b)]. Since N2 is small, the ring of com- 
ponent 2 splits into droplets, which enter the component 
1 forming small mushrooms [Fig. [2] (c)]. The droplets 
of component 2 then go towards the center and gather, 
where their complicated shapes are similar to air bubbles 
rising in water. 



Rayleigh- Taylor instability in an isotropic 
system 



Next we consider a system confined in an isotropic trap 
given by Vj = mjOJ 2 r 2 /2 with r 2 = x 2 + y 2 + z 2 . The 
initial state is the ground state for an = 200aB and 
uji = UJ2 1 in which component 2 with a spherical shape is 
surrounded by a shell of component 1 [Fig. [3] (a)]. The 
trap frequency uo\ of component 1 is then increased grad- 
ually. The outer component is pushed inward by the in- 
crease in the trap frequency and the RTI is induced at the 
spherical interface. At t ~ 36 ms, the RTI breaks the ro- 
tation symmetry and the spherical interface is modulated 
[Figs. [3] (b) and [3] (d)]. The interface is then deformed 
into a "mushroom ball" [Fig. [3] (e)]. 

The unstable modes of the interface is estimated by a 
simple analysis. We assume inviscid, incompressible, and 
irrotational fluids and component 2 of a spherical bubble 
with radius R is surrounded by component 1. The exci- 
tation frequency Q of the interfacial mode proportional 
to the spherical harmonics Y l m (0^ <j)) is given by [14[ 
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where rij is the atomic density, fj is the external force 
acting on an atom at the interface, and a is the interfacial 
tension coefficient. If Q is pure imaginary, i.e., the right- 
hand side of Eq. (jl]) is negative, the mode is dynamically 
unstable. Using the expression of a for a two-component 
BEC derived in Ref. [15| and fj = rrijUJ 2 R, we find that 
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FIG. 3: (Color online) (a)-(c) Dynamics of the cross-sectional 
density profiles d\ — \ipi(z — 0)| 2 and d 2 — \^2{z — 0)| 2 of 
components 1 and 2 and (d), (e) the isodensity surfaces of 
component 1 in an isotropic trap with frequency uji (t = 0) = 
UJ2 — 27T x 33.3 Hz. The trap frequency cji is increased such 
that uj\ is linearly increased from (c^i/c^) 2 = 1 to 3 between 
t = and t — 30 ms, and after that (cji/cc^) 2 is fixed to 3. 
The scattering length of component 1 is an = 200as and the 
numbers of atoms are Ni = N2 = 5.2 x 10 6 . The unit of the 
density is 3.0 x 10 14 cm -3 . The field of view of each panel is 
56.6 x 56.6 jLtm. 



the modes for 1 < I < 7 are unstable for the parameters 
in Fig. O which seems to be consistent with the interfa- 
cial pattern shown in Fig. [3j However, this estimation is 
only qualitative since the compressibility and the inho- 
mogeneous density distribution of a BEC are not taken 
into account. 

From Eq. (j4j), we find that the RTI is induced by an 
increase in pi or f\ , or by a decrease in p2 or ji • The den- 
sity pj depends on the interaction: an increase (decrease) 
in djj expands (contracts) component j, which decreases 
(increases) pj. The force fj acting on each component 
can be controlled, if the external trapping potential for 
each component can be controlled independently. The 
RTI can thus be induced in several ways: for example, 
(i) an increase in the scattering length of inner compo- 
nent, (ii) a decrease in the trap frequency for inner com- 
ponent, (iii) a decrease in the scattering length of outer 
component, and (iv) an increase in the trap frequency for 
outer component. The dynamics shown in Figs. [T] and [3] 
correspond to (i) and (iv), respectively. We have numer- 
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FIG. 4: (Color online) Imaginary part of the Bogoliubov 
excitation frequency, Im Q, as a function of A = (cji/cj2) 2 - 
The parameters are the same as those in Fig. [3] The modes 
for I < 8 are plotted, where / is defined in Eq. J5}. 
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FIG. 5: (Color online) Dynamics of the column density pro- 
files D\ and D2, in which an is oscillated as in Eq. © with 
ao = 80aB, A = 0.4, and cj a = 2ty x 30 Hz. Other parameters 
are the same as those in Fig. [1] The unit of the column den- 
sity is 10 12 cm -2 . The field of view of each panel is 65.4 x 65.4 
jLtm. 



C. Oscillation of the interaction 



ically confirmed that the RTI can be observed by all the 
methods (i)-(iv) for both axisymmetric oblate trap and 
isotropic trap. 

For more precise understanding of the instability, we 
perform the Bogoliubov analysis for an isotropic trap. 
We expand the GP equation (|3]) up to the first order 
of the deviation Sipj(r) from the initial stationary state 
^j(r) with spherical symmetry. The excitation mode of 
the form 

Siij = Uj (r)Yr(e, <j>) + v*(r)Yr(0, <t>) (5) 
obeys the Bogoliubov-de Gennes equation (j ^ j') 

(K jt + Vj - fij + 2g jj V 2 j + gjj^l) Uj + gjj^vj 



(6a) 



where fij is the chemical potential and 
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The stationary wave function \£j is assumed to be real 
without loss of generality. We numerically diagonalize 
Eq. ([6]) to study the stability of the system. If there is 
a complex frequency fJ, the corresponding mode grows 
exponentially and the system is dynamically unstable. 
Figure |4] shows the imaginary part of the Bogoliubov 
excitation frequency, Im as a function of (cji/c^) 2 - 
The critical value of (0J1/0J2) 2 above which Im Q rises 
increases with an increase in /, and above this critical 
value of (CJ1/CJ2) 2 , Im Q monotonically increases. At 
(0J1/UJ2) 2 = 3, which corresponds to Fig. [3j Im Q for 
the I = 5-7 modes are comparably large and these modes 
dominate the unstable dynamics. 



We demonstrate another symmetry breaking dynam- 
ics that is not due to the RTI. We return to the tight 
pancake-shaped trap used in Fig. [TJ and study the dy- 
namics for an oscillating interaction. The scattering 
length of component 1 is oscillated as 



an = ao(l + Asinuj a t), 



(8) 



where ao, A, and uo a are constants. Figure [5] shows the 
dynamics of the column density profile. By the oscil- 
lating repulsive interaction of the inner component, the 
circular interface undergoes breathing oscillation. The 
axisymmetry of the interface is then spontaneously bro- 
ken [Fig. [5] (b)], which is understood as the parametric 
amplification of an interface mode. The interface subse- 
quently exhibits various patterns [Figs. [5] (b) and[5](c)]. 
A variety of patterns can be observed depending on the 
parameters in Eq. (|5J) (data not shown). 



IV. CONCLUSIONS 

In conclusion, we have investigated the interfacial in- 
stabilities and subsequent dynamics in phase-separated 
two-component BECs. Since the initial state has ro- 
tation symmetry, the symmetry breaking nature of the 
RTI can specifically be observed in this system. We have 
demonstrated the RTI and ensuing dynamics for an ax- 
isymmetric oblate trap (Figs. Q] and [2j) and an isotropic 
trap (Fig. [3]), and the mushroom-shaped patterns are ob- 
served for both systems breaking the rotation symmetry. 
We performed the Bogoliubov analysis for the isotropic 
system and obtained the unstable spectrum (Fig. [4]). We 
also examined the dynamics for oscillating interaction 
and found that the axisymmetry is spontaneously bro- 
ken and nonlinear patterns emerge (Fig. [SJ. 
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In view of the recent development in the control of two- 
component BECs [l3|, EH, 03 ? we expect that not only the 
phenomena predicted in the present paper but also other 
theoretical predictions [l8l-[23j concerning the interfacial 
instabilities in two-component BECs will be realized in 
experiments in the near future. 
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